A New Unified Dark Fluid Model and Its Cosmic Constraint 



Lixin XCl 

Institute of Theoretical Physics, Dalian University of Technology, Dalian, 116024, Peoples Republic of China 

In this paper, we propose a new unified dark fluid (UDF) model with equation of state (EoS) 
w{a) = — a/(/3a~" + l), which includes the generalized Chaplygin gas model (gGg) as its special case, 
where a, j3 and n are three positive numbers. It is clear that this model reduces to the gCg model 
with EoS w{a) = -Bs/{Bs + il-Bs)a-^^^+°''>), when a = 1, /3 = {1-Bs)/Bs and n = 3(1 + a). By 
combination the cold dark matter and the cosmological constant, one can coin a EoS of unified dark 
fluid in the form of w{a) = —1/(1 + (1 — QA)a~^ /Qa)- With this observations, our proposed EoS 
provides a possible deviation from ACDM model when the model parameters a and n deviate from 
1 and 3 respectively. By using the currently available cosmic observations from type la supernovae 
(SN la) Union2.1, baryon acoustic oscillation (BAO) and cosmic microwave background radiation 
(CMB), we test the viability of this model and detect the possible devotion from the ACDM model. 
The results show that the new UDF model fits the cosmic observation as well as that of the ACDM 
model and no deviation is found from the ACDM model in 3a confidence level. However, our new 
UDF model can give a non-zero sound speed, as a contrast, which is zero for the ACDM model. We 
expect the large structure formation information can distinct the new UDF model from the ACDM 
model. 



I. INTRODUCTION 

The cosmic observations from type la supernovae (SN 
la) indicate that our Universe is undergoing an accel- 
erated expansion Jj, JJ. Holding the Einstein's gravity 
theory, the observed energy component T°^^ and the ge- 
ometric structure of our Universe in hand, one can deduce 
the existence of a remained un-obsered energy component 



T, 



dark 



rpobs 



(1) 



which is dubbed as a dark fluid (DF). Concerning this 
DF, in the literature, one usually separates it into two 
energy components, i.e. the named dark matter and dark 
energy for recent reviews on dark energy please see , 
or treats it as a unified dark fluid (UDF) [Tol - [l9| . 

The most popular UDF models include the so-called 
Chaplygin gas (Cg) model, its extensions the gCg model 
and the modified Chaplygin gas (niCg) model and the 
model with a constant adiabatic sound speed (CASS), 
which have been discussed extensively in the literature 
where cosmic observations from SN la, BAO, 
CMB, Gamma-ray bursts and OHD data points have 
been used to test the viability of this UDF model and 
to constrain the model parameter space. Recently, tight 
constraints to the gCg, mCg and CASS models were ob- 
tained by using the data points from SN la, BAO and 
fuU CMB information from WMAP7 ^MM- ^ was also 
pointed out that they are competitive cosmological mod- 
els to ACDM model ^MM- 

For a perfect fluid with EoS w{a) — p/ p, one can easily 
obtain the adiabatic sound speed cj,{a) — Sp/Sp — p/p in 
the form of 



9 , , , , 1 dln(l + w(a)) 



din a 



(2) 
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which characterizes the propagation of linear perturba- 
tion of UDF. When a EoS is specified, one can derive 
the corresponding adiabatic sound speed, and vice versa. 
In fact, for a constant adiabatic sound speed case, the 
model has been investigated in Ref. [llj. A general case 
was also discussed in [20], where a relation = a{—w)'^ 
was proposed. In this model, when j — 1, it reduces to 
the gCg model. As known for a gCg model, the EoS is 
given as 



w{a) = -B,/{B, + {l-Bs 



-3(H-a) 



(3) 



It would be interesting to consider the EoS of gCg from 
another points of view. Let us consider the ACDM model 
at first, one can composite the cold dark matter and the 
vacuum energy into a UDF, then the coined UDF has 
EoS 



w{a) 



f^A + (1 - ^A)a 



-3 ■ 



(4) 



Along this line, for a quintessence dark energy model 
with a constant EoS w, one can also coin a UDF with 
EoS 



w{a) = — 



fideO + (1 - f^<ieo)a3"' ' 



(5) 



Put it clearly, one can rewrite the Eqs. ^ [S]) into a 
uniformed form 



w{a) = — 



1 



/3a-" + 1 ■ 



(6) 



where /3 = (1 - Bs)/Bs or /3 = (1 - ndeo)/^deO i^deo = 
r^A) and n = 3(1 -I- a) or n = —3w {w = —1) for the gCg 
and quintessence (ACDM) model respectively. Then one 
can see that the gCg model is nothing but from a coined 
EoS from the ACDM model. So one would not be sur- 
prised that the currently available geometric observations 
from SN la, BAO, CMB cannot distinguish the ACDM 
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model from the gCg model. This situation becomes more 
serious when one uses the derived values of the CMB shift 
parameters (z*, la and R) based on ACDM model. This 
circular problem is the obvious drawback when the de- 
rived CMB shift parameters are used to test viability or 
constrain the model parameter space. 

Along this line, for generalization, we propose a new 
UDF model with the EoS given in the form of 



w{a) = -a/(^a-" + 1), 



(7) 



where a, /3 and n are three positive numbers. This EoS 
includes the gCg model as its special case when a = 1, 
/3 = (1 — Bs)/Bs and n = 3(1 + a). A novel character 
of this EoS is the fact that it gives more possibilities to 
deviate from the ACDM model in a simple way when a 
and n deviate from 1 and 3 respectively. 

In this paper, we do not want to decompose the UDF 
energy component into the so-called dark matter and 
dark energy components. As we have pointed out in our 
previous papers [Tll - fTsj that the decomposition is not 
unique. Then unphysical effect would be introduced by 
decomposition. 

Then one can use the currently available cosmic obser- 
vations to test the viability of this model and investigate 
the possible deviation from ACDM model. When one is 
doing the test, the more important is to avoid the cir- 
cular problem. To satisfy this condition, we will use the 
SN la, BAG and full information of CMB from WMAP7 
instead of the derived CMB shift parameters to test the 
viability of the new UDF model. 

This paper is structured as follows. In section |lTl the 
basic equations for UDF are shown which include the 
background and perturbation equations. In section IIII[ 
the constraint method and results are presented. A sum- 
mary is given in Section [IVI 



Then one has the Friedmann equation from the Einstein 
equations 

+ nu [(a" + /3)/(l + /3)]'"/" a-^} (10) 

where H is the Hubble parameter with its current value 
Hq = 100/ikm s^^Mpc"^, and fli {i = b,r,k) are di- 
mensionless energy parameters of baryon, radiation and 
effective curvature density respectively. Here flu — 
I — Qb — — is the dimensionless energy density for 
the UDF. In this paper, we only consider the spatially 
flat FRW universe. 

The perturbation equations for UDF should be in- 
vestigated when one considers the effects on the CMB 
anisotropic power spectrum. Here, we will not consider 
any direct interaction between UDF with any other en- 
ergy components. It means that the UDF interacts with 
the rest of matter purely through gravity. By the assump- 
tion of pure adiabatic contribution to the perturbations, 
the speed of sound for UDF is written in the form 



c2(a) = w{a) 



1 dln(l 
3' 



w{a)) 



din ( 



/3n 



(1 - a)a" -f l3 



Pin - 3a) 



(11) 



which should be in the range [0, 1] to keep the perturba- 
tion evolution stable and to avoid casual problem. How- 
ever, it is not easy to give the parameter space analyti- 
cally where is a small positive number for any value of 
scale factor a. We will use the currently available cosmic 
observations to constrain the possible parameter space. 
We will explain the technical issue in the next section. 

In the synchronous gauge, using the conservation of 
energy-momentum tensor T^f.^ — 0, one has the pertur- 
bation equations of density contrast and velocity diver- 
gence for UDF 



II. BASIC EQUATIONS OF COSMOLOGY 

For a space with uniform curvature, the line element is 
given by the Friedmann- Robertson- Walker (FRW) metric 



1 - fcr2 



(8) 

where k = 0, ±1 is the three-dimensional curvature and 
a{t) is the scale factor. With our proposed new EoS (O, 
considering the energy conservation, one has the energy 
density pu for the UDF 



pu = Puo exp 



Puo 



1 + /3 



[1 -I- w{a)] din a 



(9) 



5u = -{I + w){eu + -) - m{cl ~ w)5u (12) 



3u = -H{1 - 3cl)eu 



1 + w 



k'^Su - k^au (13) 



following the notation of Ma and Bcrtschinger j21i] . The 
shear perturbation au = is assumed and the adiabatic 
initial conditions are adopted in our calculation. 



III. CONSTRAINT METHOD AND RESULTS 
A. Method and data points 

We test the viability of this model by performing the 
observational constraints on parameter space via the 
Markov Chain Monte Carlo (MCMC) method which is 
contained in a publicl y av ailable cosmoMC package [1^ , 
including the CAMB [2j] code to calculate the theoreti- 
cal CMB power spectrum. We modified the code for the 
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UDF model with its perturbations included. The follow- 
ing 8-dimensional parameter space is adopted 

P = {cjb, Qs, T, a, /3, n, n„ log[10i%]} (14) 

where ujb = flbh^ is the physical baryon density, 05 
(multiplied by 100) is the ratio of the sound horizon and 
angular diameter distance, r is the optical depth, a, /3 
and n are three newly added model parameters related 
to UDF, Us is scalar spectral index, is the ampli- 
tude of of the initial power spectrum. Please notice that 
the current dimensionless energy density of UDF ilf/ is 
a derived parameter in a spatially flat (fc = 0) FRW 
universe. So, it is not included in the model parame- 
ter space P. The pivot scale of the initial scalar power 
spectrum kso = 0.05Mpc~^ is used. We take the fol- 
lowing priors to model parameters: G [0.005,0.1], 
Qs e [0.5,10], T e [0.01,0.8], a e [0,5], /? e [0,1], 
n e [0,5], Us € [0.5,1.5] and \og[10^° As] € [2.7,4]. In 
addition, the hard coded prior on the comic age lOGyr < 
to < 20Gyr is imposed. Also, the weak Gaussian prior 



on the physical baryon density uJb = 0.022 ± 0.002 [24| 
from big bang nucleosynthesis and new Hubble constant 
Ha = 74.2 ± 3.6kms"^Mpc"^ [H] are adopted. 

As is shown in Eq. dill) , the expression of c^, which 
contains model parameters a, /?, n and scale factor a, is 
complicated. Giving an explicit range of model parame- 
ters to keep nonnegative is really difficult. Maybe, in 
some senses, it is impossible. To circumvent the problem, 
we take the code as a black box, and hard code the con- 
dition 1 > Cg{a) > in the sampling process. It means 
that, in very sampling, the values of a, /3 and n can go to 
the next calculation stage only if they give the values of 
Cg satisfying the constrained condition for every value of 
the scale factor a. We can check the evolution of with 
respect to scale factor a once the final result is obtained. 
As is shown in the following section, please see Figure [5J 
the strategy really works. 

The distribution of parameter space is given by cal- 
culating the total likelihood C oc here the is 
given as 



X 



XCMB 



1 i 2 

Xbao + XSN- 



(15) 



The CMB data include temperature and polarization 
power spectrum from WMAP 7-year data [26]. The ad- 
ditional geometric constraint comes from standard ruler 
BAG and standard candle SN la. For BAG, the val- 
ues {rs(Z(i)/£'y(0.2), rs(2;d)/£'y(0.5)} and their inverse 



covariant matrix [23] are used. To use the BAG informa- 
tion, one needs to know the sound horizon at the redshift 
of drag epoch z^. Usually, is obtained by using the 
accurate fitting formula 28] which is valid if the matter 



scalings pb oc and pc oc are respected. Gbviously, 
it is not true in our case. So, we find z^ numerically from 
the following integration |29| 



dri'-Td 



V 



da R 



1 



(16) 



where R — ipb/Ap-^, ut is the Thomson cross-section and 
Xe{z) is the fraction of free electrons. Then the sound 
horizon is 



Tsizd) = / dr/Csil + z). 
Jo 



(17) 



where = 1/ 
substitution [2£ 



R) is the sound speed. We use the 



dz— — rrsizdj, 



' fs{zd) 



(18) 



to obtain unbiased parameter and error estimates, where 
dz = fsjzd) I Dv(z), fs is evaluated for the fiducial cos- 
mology of Ref. I27ll . and Zd is obtained by using the 
fitting formula [2^ for the fiducial cosmology. Here 
Dy{z) = [(1 + zfD\cz/H{z)Y/^ is the 'volume dis- 
tance' with the angular diameter distance Da- The 580 
Union2.1 data with systematic errors are also included 
3C]. For the detailed description of SN, please see Refs. 

sc.lMI. 



B. Fitting Results and Discussion 

We generate 8 independent chains in parallel and 
stop sampling by checking the worst e-values [the vari- 
ance(mean)/mean(variance) of 1/2 chains] i? — 1 of the 
order 0.01. The fitting results of the model parameters 
and derived parameters are shown in Table. HI where the 
mean values with Icr, 2cr and Bcr regions from the combi- 
nation WMAP+BAO+SN are hsted. Gorrespondingly, 
the contour plots are shown in Figure [T] For compari- 
son, using the same data combination, we show the cor- 
responding results for ACDM model also in Table HI For 
intimately touching to the ACDM model, the a = \ case 
is also explored for the same data combination. 
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Prameters NUDF Mean with errors NUDF(q = 1) Mean with errors ACDM Mean with errors 



e 

T 

a 
n 



log[10i°^.] 
Age/Gyr 

Ho 



n n9'?n+0-000562+0. 00112+0. 00176 r, 0998 
'^•'-'^'^'-'-0.000567-0.00111-0.00170 U.U^ZO 



1 0485+'''°°^'**^'''°°^°^"'"°'°'''''''' 

3+0.00657+0.0250+0.0436 



-0.00721-0.0239-0.0388 



-0.00253-0.00491-0.00763 

0.0903 

1 nnnc;/1R+0. 000863+0. 00229+0. 00427 
l.UUUd'lU.Q 000825-0.00143-0.00206 
r, O77+0. 0257+0. 0557+0. 0890 
"J-^' ' -0.0256-0.0475-0.0679 
o rin/lS7+0-00108+0-00'i'8''+0-0158 
O.UU'lO I -0.00487-0.00487-0.00487 
n QCQ+0. 0151+0. 0308+0. 0502 
U.»OO_0 0150-0.0290-0.0453 
q 0000+0.0346+0.0708+0.110 
O.UOZO_o 0346-0.0671-0.0998 
n QCC+0. 00164+0. 00323+0. 00488 
U.ydU_o 00165-0.00338-0.00530 
iq fjc;^ +0.112+0.221+0.330 
lO.UUI.o. 112-0. 220-0. 333 
n n/1/1+0-0017+0. 0034+0. 0053 
U.U44_0. 0016-0. 0032-0. 0049 
in C/l^+l. 191+2. 329+3. 670 
lU.Oy:'±_iig7_2. 403-3. 789 
79 /IQQ + 1.587+3. 140+4.900 
' ^■^'^'^-1.579-3.0947-4.596 



+0.000554+0.00111+0.00173 
-0.000561-0.001081-0.00160 



1 n/1Sc;+0-00251+0. 00484+0. 00736 
l.U^OO_o. 00257-0. 00503-0. 00784 
n nQnc;+0 00639+0.0249+0.0423 
U.UMUO_o. 00730-0. 0234-0. 0380 



n OQn+O-0258+0. 0544+0. 0887 
U.ZOU_o. 0261-0. 0476-0. 0681 
q r|r|^ir|+'' 000867+0. 00649+0. 0135 
O.UUy:lU_o 00410-0.00410-0.00410 
n QQ9+0. 0146+0. 0304+0. 0488 
U.yoz_o 0149-0.0285-0.0400 
q nsi R+0-0342+0. 0699+0. 105 
O.UOiD_o. 0346-0. 0652-0. 0949 



0.956 



+0.00168+0.00320+0.00477 



0.00168-0.00331-0.00512 
iq (.0/1+0.112+0.218+0.335 
10.00^_0iii_o. 217-0. 334 

n n/i/i+o ooi''+o oo33+o.oo5i 

U.Uy:4_0. 0017-0. 0032-0. 0048 
in CIR+1. 192+2. 340+3. 516 
lU.010_iigo-2. 299-3. 605 
79 91 1+1.556+3.118+4.898 
I ^■^11^1,546-2.983-4.461 



n n99R+0-000514+0. 00100+0. 00150 
U.UZZU_0 000509-0.00101-0.00145 
n 1 in+000367+0. 00736+0. 0113 
U.llU_o. 00368-0. 00711-0. 0110 
1 n'?QS+0-00253+0. 00493+0. 00755 
l.UOaO_o. 00252-0. 00506-0. 00755 
n nQ97+0-0064''+0-0246+0.0452 
U.uyz ( -0.00719-0.0235-0.0379 



n Q71 +0.0119+0.0236+0.0370 
U.»( 1-0.0117-0.0232-0.0342 
q r,ocq+0-0344+0. 0703+0. 106 
O.UOOO_o. 0349-0. 0664-0. 0978 
n 7/19+0.0162+0.0296+0.0440 
<1^_0. 0160-0. 0333-0. 0512 
iq 71 o+0. 106+0. 208+0. 311 
10. / lo_o 105-0.201-0.312 
n oco+0. 0160+0. 0333+0. 0512 
U.ZOO_o 0162-0.0296-0.0439 
10 g91+l l*9+2-338+3.603 



71.629 



-1.192-2.371-3.578 
+ 1.446+2.783+4.381 
-1.435-2.845-4.0537 



TABLE I. The mean values of model parameters with 1 — 3a errors from the combination WMAP-f BAO+SN. 



r 



From the Table |T] and Figure [1] one can clearly see 
that a tight constraint to the model parameters a, /3 and 
n is obtained. The results indicate that the proposed 
UDF model do not deviate from ACDM model at the 
background level distinctly. The values of a and /3 are 
very close to the ACDM model limit 1 and 3 respectively. 

Using the mean values of model parameters, we plot 
the evolution of EoS w{a) and the speed of sound c^(a) 
for UDF with respect to the scale factor a in Figure [2] 
From the left panel of Figure [21 the UDF behaves like 
cold dark matte at early epoch and dark energy at late 
epoch. The right panel of Figure [D shows that the value 
of the sound speed Cj(a) of UDF is small positive num- 
ber and varies with scale factor. The small values of c^(a) 



make it possible to form large scale structures in our uni- 
verse. Also, the positivity of c^(a) is really guaranteed 
by the 'filter' in sampling. One can read the fact that the 
Cg(a) for UDF evolves with respect to the scale factor a 
but that is constant and zero for the ACDM model from 
the Figure [2] It implies that the observational effects 
that can be used to discriminate the UDF model from 
ACDM model are due to the sound speed and EoS. The 
large structure formation information and late time in- 
tegrated Sachs- Wolfe (ISW) effect can be combined with 
background evolution data sets to discriminate the mod- 
els. However, to use the observed large structure infor- 
mation, one has to investigate the large scale structure 
formation process in UDF model which is out the range 
of the current paper. 



r 



For comparison, we show the CMB Cf^ power spec- 
trum v.s. the multiple momentum I in Figure [3] where 
the mean values of the relevant cosmological and model 
parameters are adopted from the Table U . It is clearly 
that the curves of the power spectra of UDF match up 
to that of ACDM model. The tiny difference comes from 
the hight of the fist peak and the trough at the large scale 
{I < 100). The difference can be understood easily, if one 
writes down the correspondence /3 = (1 — ildeo)/^deO- 



Large values of /3 mean large ratio of the matter and dark 
energy. Large values of ftmo = 1 — ^deo make the equality 
time of matter and radiation earlier if they evolve with 
the same scaling law, say oc a~^, then the first peak of 
the CMB power spectrum is shifted to a higher position. 
At the large scale {I < 100), the difference is due to the 
tiny deviation from the scaling law oc a~^, i.e. through 
the integrated Sachs- Wolfe (ISW) effect. 



IV. SUMMARY 



combination of the cold dark matter and the cosmological 



In summary, we proposed a new unified dark fluid 
model w{a) = —a/{(3a~"' + 1) which is inspired by a 
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FIG. 1. The 2-dimensional contours and 1-dimensional probability distribution of model parameters with 1 — 3cr errors. Bottom 
panel show the case where the mode parameter a is fixed to 1. 
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FIG. 3. The CMB C/ ^ power spectrum v.s. multiple moment I. The relevant model parameters are adopted to the mean 
values as listed in Table |I] for NUDF and ACDM model, where the black dots with error bars denote the observed data with 
their corresponding uncertainties from WMAP 7-year results, the blue dashed line is for NUDF with a free, and the green 
dashed line for NUDF with a — 1, the red solid line is for ACDM model with mean values for the same data points combination. 
And the red doted line is for ACDM model with mean values taken from i32^ with WMAF+BAO+iJo constraint results. 



model constant. In this simple form, one can see that the 
EoS of the gCg model is nothing but a coined one from 
that of ACDM model. And the gCg model is a special 
case of our proposed one. The new unified model gives 
more possibilities to deviated from the ACDM model in 
the unified dark fluid case when the model parameters a 



and n take different values a and 3 respectively. To test 
the viability of this new UDF model and detect the pos- 
sible deviation from ACDM model, we perform a global 
fitting by using the MCMC method with a combination 
of the full CMB, BAG and SN la data points. The re- 
sults show that the new UDF model fits the observational 
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data sets as well as that of ACDM model. And, one does 
not see any deviation from ACDM model in Sc confidence 
level under the currently available data sets. But I would 
like to point out that the new proposed UDF model can 
give a non-zero sound speed. As a contrast, the coined 
EoS from ACDM model has zero sound speed. Base on 
this point, we expect large scale structure formation in- 
formation can distinct the new UDF model from ACDM 



model. 
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